Konstanz University • Theory Group 
( Gravity | Quantum Theory | Optics ) 



Preprint KONS-RGKU-95-01 - \quant-ph/950301 



in 

On 



Relation between energy shifts and relaxation rates 
for a small system coupled to a reservoir 

■ Jiirgen AudretschQ, Rainer MiillerQ and Markus Holzmann 

Fakultdt fiir Physik der Universitdt Konstanz 
^ ■ Postfach 5560 M 674, D-78434 Konstanz, Germany 

O ' 



in 

On 



Abstract 

For a small system the coupling to a reservoir causes energy shifts as well as transitions 
between the system's energy levels. We show for a general stationary situation that the energy 
shifts can essentially be reduced to the relaxation rates. The effects of reservoir fluctuations 
and self reaction are treated separately. We apply the results to a two-level atom coupled to 
a reservoir which may be the vacuum of a radiation field. 
^ ' PACS numbers: 32.80-t; 42.50-p. 

cr 

> 

1 Introduction 

For a small quantum system which is coupled to a reservoir, energy shifts due to the coupling will occur 
as well as transitions between the system's energy levels. Our aim in the present letter is to derive 
relations between these relaxation rates and energy shifts which are valid under very general conditions. 
We consider a system which moves in a stationary way on a possibly accelerated trajectory in a general 
spacetime (which may be curved or possess nontrivial boundaries). The stationarity of the situation 
demands that the system moves along the orbits of a timelike Killing vector field. Stationarity is also 
assumed for the reservoir. The case of an extended system at rest is included. Typical realizations of 
such a situation are an atom or an elementary particle (system) coupled to a quantized radiation field in 
the vacuum or a many-photon state (reservoir), whereby the atom may be accelerated. In this case the 
energy shift is the Lamb shift. 

For the description of the system-reservoir interaction we will generalize the formalism which has 
been established by Dalibard, Dupont-Roc and Cohen- Tannoudji (DDC) Q and extended in |^|, [| to 
the situation just described. It leads to a clear and physically appealing interpretation of the processes 
in terms of the underlying physical mechanisms since it allows the separate discussion of the effects of 
reservoir fluctuations and self reaction (or radiation reaction). 
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2 Reservoir fluctuations and self reaction 



First we provide the necessary general relations. We consider a spacetime which is covered by a coordinate 
system x = (t, x). The time coordinate t is assumed to be the "natural" time for the description of the 
reservoir. The system is described by a stationary (accelerated) trajectory x(t) — (t(r), x(r)). This may 
be realized for example by a pointlike object moving on this trajectory. The respective proper time 
variable t is used for the parametrization of the trajectory. 

The time evolution of system and reservoir has to be specified in terms of a single time variable for 
which we choose the system time r. The Hamiltonian which governs the dynamics with respect to r is 
given by 

H = Hs(r)+H a (t)^+V(x(r)) (1) 

where H$ is the free Hamiltonian of the system. Hn(t) is the free reservoir Hamiltonian with respect 
to the time variable t. The factor dt/dr in ([!]) is due to the corresponding change of variables. V{x{t)) 
represents the coupling between system and reservoir. It is effective only on the trajectory x(t) of the 
system. It is assumed to have the general linear form 

V(x(T)) = -gJ2Ri(r)Si(x(r)) (2) 

i 

with g being a coupling constant. Ri and Si are hermitean reservoir and system operators. For an 
extended system at rest, we can set r = t and replace x(t) by t here and in the equations below. 

We can now write down the Heisenberg equations of motion for observables of the system and the 
reservoir. For our purposes, the following will be important: It is possible in the solution of the Heisenberg 
equation for the reservoir variable Ri to distinguish on one hand the part which is present even in the 
absence of the coupling. It is independent of the system and is called the free part R{ of Ri . The remaining 
contribution is caused by the presence of the system and contains the coupling constant g. It is called 
the source part Rf of R,^. 

R i (T)=R{(r) + Rl(T). (3) 

We consider the rate of change of an arbitrary system variable G. Because of the coupling (^|), reservoir 
operators appear in the Heisenberg equations of G. According to (|^), they can be divided into their free 
and source parts. The rate of change of G due to the coupling can therefore be split into two contributions 
which correspond to two different physical mechanisms: (i) the change in G produced by the fluctuations 
of the reservoir which are present even in the absence of the system - this portion is related to the 
free part of the reservoir and is called the contribution of reservoir fluctuations to dG/dr - and (ii) the 
change in G due to the interaction with the excitations of the reservoir which are caused by the system 
itself. This is the contribution of self reaction or radiation reaction and is connected with the source part 
of the field. Following DDC Q, we adopt a symmetric ordering between system and reservoir operators. 

In a perturbative approach, we take into account only terms up to second order in g. Since we are 
interested only in the dynamics of the system, we average over the reservoir degrees of freedom. We 
assume that the density matrix is factorized into a system and a reservoir part at the initial time t = 0: 
p(0) = ps(0) Pr(0)- We select one specific system state \a) and take the expectation value with respect 
to that state. 

Proceeding essentially as in [Q, we find for the contribution of reservoir fluctuations 
(^W) =i {[ H r f / (r)Mr)])-^^{[Yi(r),[S{(r),G(r)]] + [s{(r), K(r), G(r)]] ) (4) 

whereas the contribution of self reaction is (curly brackets denote the anticommutator): 

dG \ - 2 
dr 



(r)\ =i{[Htt f {r),G{r)])- 9 -Y,{{Hr), [S{(t),G(t)]} ~ {^/(r), [Z^r), G(r)]} ) . (5) 

' i 

In (|J) and (||), brackets (...) denote averaging over the reservoir and taking the expectation value in the 
system state \a). Furthermore we have introduced the effective Hamiltonians 

K"<r) % E [W' . : = -% E {Ur\ S{(r)} (6) 
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ith 



Y % (r):=Y, dr'C^(x(r),x(r'))Sf(r'), Z<(r) := £ / dr' x$(x(t), x(t'))SJ(t') (7) 
. Jo . Jo 

Since we are interested in time scales which are large compared with the correlation time of the reservoir, 
we have extended the range of integration in (Q) to infinity. We have also introduced the symmetric 
correlation function and the linear susceptibility xfi °f the reservoir: 



X K(x(t),x(t>)) :=±Tr R (p„(0) [r{ (x(t)), R? (x(r')) 



(8) 
(9) 



Because of the stationarity, C^j and Xu are only functions of the time difference t — r' . 



3 Energy shifts and relaxation rates 

The relaxation rates of the system's energy in the state \a) can be obtained from (0) and (0) with the 

f f l-i f 

choice G = Hs- We replace [S( , H J S ] in second order by i-j-S( and find for the contributions of vacuum 
fluctuations and self reaction to the system's relaxation rate 



« ( 2 (T)) z(/))^4.( T y) ( io) 



(11) 



with the symmetric correlation function and the linear susceptibility of the system 

CS(t,t>) := \{a\ {s{(r), s/(r')} |a> = £ Re (<a|S/(0)|6)<6|S/(0)|a>e*-»^ '>) , (12) 

6 

Xg(r,r') = ~ (a| [fif/ (-r) , 5/ (r')] |o> = z^Im (<a|S/ (0)|6>(6|5/(0)|a)e^ ^ '>) . (13) 

b 

To evaluate (jlj) and (|l^) we have used the stationarity of the situation which allowed us to introduce a 
complete set of stationary system states \b) (eigenstates of Hg) with energies u>b and uj a b = to a — u>i,. The 
formulas © and (£|) are the generalizations of the corresponding equations of DDC |jj to the situation 
considered here. (Note the slight differences in the definition of the correlation functions). 

Beneath that, the system-reservoir coupling leads to a shift of the system's energy levels. The second 
order radiative shift of state \a) is given by the expectation value of the effective Hamiltonians @ in 
that level. Again, the total shift can be split into the contributions of vacuum fluctuations and radiation 
reaction: 

poo 

{8E a ) rf = - ig 2 Y; dr' C^{x{r),x{r'))xUry) (14) 
i,j ■>'■ 

pOC 

(SE a ) sr = -i 9 ^ dT' X K(x(T) 7 x(Tl))CS(T,Tl) (15) 

i,j J ° 

These formulas too are generalizations of those in Qj. Because {xfj) is symmetric (antisymmetric) 
in i and j, terms with uj a b = do not contribute to the relaxation rates and energy shifts. The proof, 
which wc omit, is based on the stationarity. Without restriction, we can therefore exclude these terms 
below from all b summations. 
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4 Relation between energy shifts and relaxation rates 

We now come to the main point of this letter. We will prove to second order in g for a system in the 
state | a) quite general relations between the corresponding relaxation rates ( flO| ) and (|ll]) and the energy 
shifts (|l4|) and (|l5|). We start with the contribution of reservoir fluctuations. Using the explicit formula 
( |l3| ) for the system's linear susceptibility, the energy shift ([l4]) can be written 

(SE a ) rf = -g 2 ^ f°° dr'C^(x(r),x(r'))^Re ((a\S{(0)\b)(b\Sj(0)\a)^-e^ T - T 'A (16) 

i,j,b ^ a ' 

Now we can apply the Kramers-Kronig relation |^| 

Re(fM) = - / du/lm(jV))- (17) 



to the function f(oj ab ) = (a\S{ (0)|&)(&|S/ (0)|a) u> ab e ttAjab ^ T T > which is analytic for u Q {, ^ and vanishes 
for uj a b — > ioo. We find 

2 r + oo poo 

(SE a ) rf = <W dr'C§(x(r),x(r')) (18) 



7T 



x-^Im f(a|S/(0)|6)(6|5/(0)|a)l e -'( T -')) - 



V 



dr \ 1 J lu' J uj' — uj a b 

Introducing the quantity r r a l(ui') by the definition 



ri(c'):=-2 5 2 ^ / dr'C^(x(r),a ; (r'))-Im(( a |S/(0)|6)(6|5/(0)|a}e-'(— '>) (19) 
we obtain the desired relation 



TJ 2tt ^ 7_ oc uj' uj' - uj ab 

Eq. ([20|) connects energy shift and relaxation rate since the latter can be expressed easily in terms of 

dH s ' 
dr 



T rf. 



= J2KlM, (21) 

r f b 

where we have used Eqs. ( px| ) and (|l^). Following essentially the same procedure, an equation analogous 
to ( pfl| ) can be derived for the contribution of self reaction to the energy shift 



1 r+oo 

<'*>-- s?L 



i/SM * , (22) 



where 

poo J . 

rr b (^):-2z 5 2 ^ y dr' X ^(x(r),x(T'))^Be[(a\S{(0)\b)(b\S^O)\a)e^ ^-^) (23) 
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and Eq. (|2l] ) holds correspondingly (replace r/ by sr). Note that r^ sr do not depend on r because of 
the stationarity of the physical situation. 

We have studied above a very general stationary situation: arbitrary linear system-reservoir coupling, 
arbitrary stationary motion in a flat or curved spacetime, nontrivial boundary conditions allowed, arbi- 
trary state of the reservoir. We have shown that the determination of energy shift and relaxation rate 
can be reduced in a unified way directly to the calculation of the coefficients T r JJ sr . They turn out to be 
the fundamental underlying quantities. The calculation of the energy shift has thereby been simplified 
as compared to ( |l4| ) and (|l5|). Equations (|2(]), ( ^l| ) and ( p2| ) show that the T r JJ sr refer to particular 
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transitions \a) — * |6). This becomes even more evident when the system has only a few energy states 
as will be demonstrated in the next section. In addition we see that the mutual dependence between 
energy shifts and relaxation rates holds for the reservoir fluctuation terms and for the self reaction terms 
separately. This supports the view that the distinction of these two mechanisms is physically reasonable. 
Finally we mention that for concrete physical situations it may be necessary to introduce a frequency 
cutoff in order to regularize (5E a ) r f / sr . As compared with ( |l4| ) and (|l5|), the new expressions @ and 
(^H) are directly prepared for this. 

5 Application: Two-level atom 

To demonstrate the usefulness of the proposed scheme we consider a two-level system moving on a 
stationary trajectory and derive several statements which are valid under very general conditions. This 
illustrates at the same time the physical concepts. We assume two stationary states |+) and |— ) with 
energies ±|wo- The atomic Hamiltonian can be written with the help of the pseudospin operator S3 = 
|l+)(+|-l|-)Has 

H s = oj S 3 {t) (24) 
The coupling to the reservoir is linear and connects only different levels of the system: 

V = - 9 S 2 (t)R(x(t)) (25) 
with S 2 = |(5 + - S-) and S± = |±)(T|- We rewrite T r J b /sr from Eqs. @ and @ as 

Kl /Sr M = ~2Lu ab \{a\S((0)\b)\ 2 Y f/sr M (26) 
where we have simply \(a\S[ (0)|6)| 2 = j for a ^ b. Here, 

Y f M = g 2 / dr'C R (x(T), x(t')) cos uj ab {T - /), (27) 
Jo 

/•oo 

Y r (u ab ) = ig 2 / dr'x R (x(T), x(t')) smuj ab (T - r 1 ), (28) 







and u) a b = iu)o for a > b (a < b). 

Using the general expression (|2l]) for the relaxation rate together with (p6|), we can write 



b 

Since j r f is symmetric in ui ab , it follows that 



(^r) / = - 2 T,^\{a\Si{Q)\b)\ 2 l rf M 



(29) 



/^£)\ = -2a, o7 r/ M (E K«l4(0)|6)| 2 + E IH^(0)|fc)| 2 J (30) 

\ aT I rf \ b<a b>a ) 

Analogously, we find for the contribution of self reaction to the relaxation rate 

/^£\ = -2 W o7 sr (^o) (E K«l^(0)|6)| 2 ~ E K«l4(0)|6>| 2 ) • (31) 

\ aT ' sr \ b<a b>a ) 

From these general relations, it is possible to find expressions for the atom's evolution into equilibrium 
as well as for the Einstein coefficients Af and corresponding to upwards and downwards transitions. 
We consider the total transition rate (^) tot = (^§r) rf + (^§r) rf - Usin g © and @> [t can be 
simplified by noting 

E l(a|^(0)|6>| 2 ± E \(*\Sl(0)\b)\ 2 = { l^ (alHqla) (32) 

b>a b<a L 2uJ ° 
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Here we have replaced Hg by the total system Hamiltonian, which is justified in order g 2 . We thus obtain 
a differential equation for the evolution of the mean atomic excitation energy 

(^( T )) t = -7^o7 sr M - 7 r/ M (H s (r)) (33) 

which has the solution 

"I 1 T j SV / 1 SV \ 

(Hs(t)) = --LJo + ~ f 7 + ((Hs(0)) + 2^^jJ (34) 

A number of interesting points can be inferred from (HJ) : First the relaxation into equilibrium is always 
exponential, without oscillations, no matter what the the trajectory of the atom or the state of the 
reservoir is. Its rate is determined by the contribution of reservoir fluctuations r y r ^ alone. Furthermore, 
the first two terms of ( |34| ) show that the equilibrium excitation above the ground state — hu>o is governed 
by the relative magnitude of "f r ^ and 7 s1 *. 

The Einstein coefficients Af and Aj_ can be identified by comparison of ([mJ) with an appropriate rate 
equation (Eq. (65) of Q). One finds 

= \l rS \l S \ Al = l -Y S + \l sr - (35) 

The important consequence is that r y r ^ and 7 s1 * (and therefore T r J b and T s J b ) can separately be determined 
as functions of the Einstein coefficients. In principle they are therefore measurable quantities. Beyond 



that one can read off from (35) the physical origin of the Unruh effect |6[. Taking as reservoir the 
vacuum, the spontaneous excitation of an accelerated two-level atom results from an imbalance between 
the contributions of vacuum fluctuations and radiation reaction as has been discussed in detail in Q . 

The radiative energy shift, which in this case is called the Lamb shift of the two-level atom can now be 
obtained by application of the formulas (BO) and (|2^) . First we notice that generally the contribution of 



self reaction does not contribute to the relative shift of the two levels. From (22) it follows after a short 
calculation that 

A sr = SE s + r - SE s _ r = 0. (36) 
On the other hand, the contribution of reservoir fluctuations is 

A tot = A r/ = SE r + f - SE r J 

= -i- [°°dw'Y f (u') (-^- 7^—], (37) 

2irJ 1 v ; v^' + ^o uZ-uoJ 

where the symmetry properties of 7 r ^(w') have been used. It is determined by the fluctuations of the 
reservoir in accordance with the heuristic picture of Welton . In the derivation of Eq. (|37]) , no other 
properties of the reservoir and no further assumptions about boundaries or the atomic trajectory were 
used. Their influences are contained in the form of 7^. Note that the energy shift can be determined 
once the functional form of the Einstein coefficients (|35|) is known. This is so much the more remarkable 
since, in contrast to the usual expression for the Lamb shift, they can be calculated easily for a given 
system, for example from Fermi's golden rule. 
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